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Abstract. We introduce a coupled pair of evolution equations for the uninte¬ 
grated gluon distribution and the sea quark distribution which incorporate both 
the resummed leading ln(l/x) BFKL contributions and the resummed leading 
ln((5^) DGLAP contributions. We solve these unified equations in the perturba¬ 
tive QCD domain. With only two physically motivated parameters we obtain an 
excellent description of the HERA F 2 data. 


1. INTRODUCTION 

One of the striking features of the measurements of deep inelastic scattering 
at HERA is the strong rise of the proton structure function F 2 as x decreases 
from 10“^ to below 10““^ [1]. At hrst sight it appeared that the rise was due to 
the (BFKL) resummation of leading ln(l/x) contributions. In this approach 
the basic dynamical quantity at small x is the gluon distribution f{x,k^) 
unintegrated over its transverse momentum kp. Observables are computed in 
terms of / via the kp factorization theorem. For example 

F 2 = F^^ ® f where F]^ = X (1) 

<? 

and where Gi denotes a convolution in transverse {kp), as well as longitudinal, 
momentum, see Fig. 1(a). Fi^^ is the off-shell gluon structure function which 
at lowest order is given by the “quark box and crossed-box” contributions, 
U9 QQ- On the other hand the conventional DGLAP approach is based on 
collinear factorization 

dF 

xa/n(QVA2) = X e? Pqg®g + E + f) (2) 




FIGURE 1. (a) The diagrammatic representation of the kr factorization formula 

F 2 = (g) /. At lowest order F^^ is given by the quark box shown (together with 

the crossed box) summed over all types of quark, (b) The different regions of integration 
used to evaluate the sea quark distribution. In regions II and III the sea is driven by the 
gluon, Sq = Bq® /. 


where ® is simply a convolution over longitudinal momentum. DGLAP evolu¬ 
tion effectively sums up the leading InQ^ contributions, and is able to describe 
the F 2 data at the smallest x values observed (even for ~ IGeV^) with 
an appropriate choice of input parton distributions. These non-perturbative 
input shapes in x mean that there is more freedom in the pure DGLAP de¬ 
scription than in the BFKL approach. In the latter approach, the shape in 
X emerges, in principle, from the perturbative ln(l/a;) resummation. In prac¬ 
tice it is not so clear cut since non-perturbative, and also subleading ln(l/a;), 
effects can modify the BFKL prediction. Here we should attempt to hnd a 
unified description of F 2 which incorporates both (DGLAP and BFKL) of 
these perturbative effects. 

2. UNIFIED BFKL/DGLAP FORMALISM 

We argue that the unintegrated gluon distribution f{x, k^) and the kr fac¬ 
torization theorem provides the natural framework for describing observables 
at small x. To determine / we arrange the BFKL equation so that we only 
need to solve it in the perturbative domain k‘^ > k^ [2], We also include the 
residual DGLAP contributions. To be precise we have 
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where — 1/;^ is taken from DGLAP because it is already included in BFKL. The 
input term comes from two sources: the parts of BFKL and DGLAP 

terms. We specify the input in terms of a simple two parameter form 5 f(a;, k^) = 
N{1 — xY. In addition to restricting the solution of the BFKL equation to 
the perturbative region k'^ > k^ and to including the DGLAP terms, we have 
also introduced a 9 function which imposes the constraint k’^ < k^/z on the 
real gluon emissions. The origin of this constraint^ is the requirement that the 
virtuality of the exchanged gluon is dominated by its transverse momentum 
l/c'P ~ k'Y We take a running coupling asik^)) which is supported by the 
results of the next-to-leading order ln(l/a;) analyses of Fadin, Lipatov, Gamici 
and Giafaloni. The hnal term in (3) depends on the quark singlet momentum 
distribution S. At small x the sea quark components Sq of S dominate. They 
are driven by the gluon via the g ^ qq transitions, that is = Bq® f 
where at lowest order Bq is the box (and crossed box) contribution indicated 
in Fig. 1(a). Besides the 2 ; and k"^ integrations symbolically denoted by 0 
the box contribution implicitly includes an integration over the transverse 
momentum kt of the exchanged quark. The evolution equation for S may be 
written in the form 

S = Bq{k‘^ = zg{z,kl) + f + ^ ( 4 ) 
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where the hrst three terms on the right hand side are the “5^ 0 /” contribu¬ 
tions coming from the regions I, II, III of the k^ and integrations that are 
shown in Fig. 1(b). First, in the non-perturbative domain, region I, the u, d, s 
sea quark contribution is parametrized in the form = Cp a;“°'°®(l — x)® 
consistent with soft pomeron and counting rule expectations, where Cp is in¬ 
dependent of Q^. The constant Cp is hxed in terms of the two parameters, 
N and (3, by the momentum sum rule. In region II we apply the strong kx 
ordering approximation with Bq ^ Bq{k^ = 0) so that the k^ integration can 
be carried out to give a contribution proportional to g{x/z, k^). Finally in re¬ 
gion III we evaluate the full box contribution; this gives the main contribution 
and is responsible for the rise of F2 with decreasing x. The last two terms in 
(4) give the sea —sea evolution contribution, and the valence contribution 
V{x,QY which is taken directly from a recent parton set. The charm quark 

A more general treatment of the gluon ladder which incorporates both the BFKL equa¬ 
tion and DGLAP evolution is given by the CCFM equation [3], which is based on angular 
ordering of the gluon emissions. The angular ordered and kinematic constraints lead to sim¬ 
ilar subleading ln(l/a;) effects, but the kinematic constraint overrides the angular ordered 
constraint, except when in the large x domain [4]. 




FIGURE 2. The two-parameter description of some of the F 2 data at small x using 
f{x, k^) evaluated with (continuous curves) and without (dashed curves) the kinematic 
constraint < k\lz. 


component of the sea is given totally by perturbative QCD, since for k'^ < 
the box B^k"^ = 0 ) is hnite as 0 due to rric 7 ^ 0 . 

3. DESCRIPTION OF F 2 DATA AND DISCUSSION 

We solve the coupled integral equations (3) and (4) for the gluon / and 
the quark singlet S in the perturbative domain, /c|n > /cg. The only input is 
the gluon g^x^k^). We take /cg = IGeV^. We determine the values of the 
two input parameters by htting to the available data [1] for F 2 with x < 0.05 
and > 1.5GeV^. The continuous curves in Fig. 2 show the description of 
a sample of the data. Overall the fit is excellent; at least as good as that 
achieved in the recent global analyses. When we repeat the analysis with 
the kinematic constraint omitted we see that the description (given by the 
dashed curves in Fig. 2) is not so good and, moreover, the extrapolation of 
the gluon to a; ~ 0.4 no longer describes the WA70 prompt photon data. How 
important are the ln(l/a;) effects? First we replace the BFKL kernel in (3) 
by the standard DGLAP splitting function, Pgg. We hnd that the gluon is 
not changed by much in the HERA domain x > 10“^, compare the dashed 
and dotted curves in Fig. 3. This effect is well known. In the power series 
expansion of the gluon anomalous dimension in as/uj the coefficients of the 
2nd, 3rd and 5th terms are zero, where u is the moment variable. On the other 
hand, when we use pure DGLAP evolution for the quark singlet, as well as 
the gluon, the difference is pronounced; compare the dashed and dot-dashed 





FIGURE 3. The continuous and dashed curves come from solving (3) and (4) with and 
without the kinematic constraint. The dotted and dot-dashed curves are obtained using 
DGLAP in the gluon sector and in both the gluon and quark sectors respectively. The same 
input is used for all curves. 

curves in Fig. 3. 

We would like to conclude that we have a theoretically well-grounded and 
consistent formalism which, with the minimum of non-perturbative input, is 
able to give a good perturbative description of the observed structure of F 2 . 
Moreover the BFKL/DGLAP components of F 2 are decided by dynamics. In 
this way we have made a determination of the universal gluon distribution 
f{x, k"^) which can be used, via the kr factorization theorem, to predict the 
behaviour of other observables at small x. The predictions for F 2 (charm) and 
Fl can be found in [2]. 
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